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I. INTRODUCTION 


In recent years the great expansion in the development 
of semiconducting materials with acceptable thermoelectric 
characteristics has made them attractive for, in special 
applications, direct energy conversion and heat pumping. 

Thermoelectric power sources and heat pumping devices 
can offer significant advantages with respect to long life, 
years of unattended operation, freedom from environmental 
restraints, and unusually favorable specific power ratios. 

It is true that even now thermoelectric cooling and 
power sources are not serious competitors of powerful heat 
pumps and conventional power apparatus, but they are 
irreplaceable in those cases where the overriding require- 
ments are small size, simplicity of construction, quiet 
operation, weight, etc. 

Various microrefrigerators, used in radio engineering, 
medicine and industrial design, have received proper 
recognition. Fueled portable thermoelectric generators, 
especially for military applications, are already used. 

In the present work an attempt is made to design, at 
least theoretically, a thermoelectric air cooling device 
servicing the passenger compartment of a small automotive 
vehicle, in conjunction with a thermoelectric generator 
using the waste heat of the exhaust gasses of the engine 


to provide a portion of the required electrical power. 





The feasibility of the replacement of the conventional 
alternator by a thermoelectric generator is also examined. 

The introduction of the thermoelectric generator will 
improve the coefficient of performance of the thermoelectric 
cooler, and will also increase the initial installation cost. 

The design was developed from the macroscopic point of 
view, considering the present state of the art of thermo- 
electric materials. 

The total work is divided into five sections: 

1) The heat flow analysis for the semiconductor material 
under a temperature gradient and internally generated heat 
due to the flow of electric current. The results of this 
analysis are almost common to the thermoelectric cooler and 
generator. In this section a differential equation, which 
describes the heat balance in the thermoelement has been 
derived. In this equation, the cross section of the thermo- 
element is treated as a function of the length and the 
properties of the material, as thermal conductivity, Thomson 
coefficient, etc., are treated as functions of temperature. 

2) The analysis and design of the thermoelectric cooler. 
A detailed analysis of the T.E. cooler is made for both modes 
of operation, maximum heat pumping rate and maximum coeffi- 
cient of performance. The analysis is closed with comparison 
between the two modes and the selection of max COP as the 
most suitable for the present application. 

The analysis is followed by an indicated design of a 


cooler with a capacity of 7000 BTU/hour which, from a quick 





estimation, was found to be satisfactory for the selected 
environmental conditions. 

3) The analysis and design of the thermoelectric 
generator. The analysis is conducted from the standpoint 
of maximum efficiency mode. It is followed by selection of 
the proper thermoelectric materials for the assumed conditions. 
Computer techniques are employed for estimation of the 
memperature relations of the thermoelectric material proper— 
ties. The analysis is closed with the estimation of the 
available heat power from the exhaust gasses. The attempted 
design concerns a segmented type thermoelectric generator 
which is carried out exclusively by a computer program. 

4) Conclusions and recommendations for further study. 

5) Appendices 

Generally, the proposed configuration would be as it is 


shown in Figure l. 
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Engine 







T.E. Generator 


FIGURE I 


DIAGRAMMATIC ARRANGEMENT OF THE PROPOSED SYSTEM 


ll 





II. STEADY STATE HEAT FLOW IN THE SEMICONDUCTOR ELEMENT 


The derivation of the heat flow balance equation through 
the thermoelement is presented when the cross section of the 
element is a function of the length and the characteristics 
of the material, electrical resistivity and thermal conduc- 
tivity are functions of temperature. 

Under the proposed geometry of the semiconductor, the 


one dimension heat flow is evident. 





FIGURE 2 


PROPOSED GEOMETRY OF THERMOELECTRIC SEGMENT 
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The semiconductor is considered to be adiabatically 
insulated on all sides except side ABCD and EFGH. 


The temperature of the solid is assumed to be Ty at the 


front side (ABCD) and T., at the rear side (EFGH), when T 


Cc 


The characteristics of the material 


H 


is greater than To: 


are considered to be variable and more specifically to be 


functions of temperature which can be expressed as: 


ie conductivity ,(1) = a, + at + a5T 
es ee _ 2 
2) Resistivity ¢ (7) = bo ca b,T “F boT 
3) Seebeck coefficient q(T) = co + c,T + cot? 
4) Thomson coefficient 6(T) =T oat) = c,T + 205T° 


Considering the part of the semiconductor which is 
confined between the cylindrical surfaces of radii r and 
rtdr, let q(r) and q(r+dr) be the rates of heat flowing in 
and out of the differential element. 

If a current I is flowing through the semiconductor 
Joule heat will be generated due to the electrical resistance 
of the semiconductor and heat will be released or absorbed — 
depending on the direction of the flowing current — due to 
the Thomson effect. The differential of the Joule heat 


generation rate is given as: 


: (CT) | aie 


> 


uns) 





where A(r) is the surface of the differential solid element 
at radius r. 
Also the differential form of the Thomson effect is 


given as: 


7. = le(T)dT (1-2) 


Consider now, for balance the rate of heat flow out of 
the differential solid element must be equal to the rate of 
heat flow in, plus the rate of heat generation in the element. 


These can be expressed as: 


q (ede) = Oa a (1-3) 


But q(r+dr) can be expanded in a series of the form: 
qferde) = 96) + ao 
Or 


If the first two terms will be considered, then equation 


(1-3) can be written as: 


Wd + “ 
sp de = 4,4 Cas 


T 


Substituting the equivalent relations for qs and Qn equation 


(1-5) becomes 


dVr) _ eC) + ar , 
= > ie I Os (1-6) 


14 





The conduction heat — rate of flow — q(r) is given by: 
q(v)= - yA (1-7) 
Tr 


By differentiation in respect to r equation (1-7) becomes: 


die. 4 cn ale Arye ~y piel Uo eT \oAns a) 


By substituting of equation (1-8), equation (1-6) becomes: 


\nAwet 


pa 


4d A) | era d[A@] dT nites CAMEL : 


For the geometric configuration of Figure 2 the cross 


section area A(r) is a linear function of radius r: 


Ar) = eter (1-10) 
A(r) = me (1-10a) 


where m =PT ‘ 


Using the above expression for the cross section, 


equation (1-9) becomes: 


xm mod T,sbo] met]. [aan s Teens + EOP. 0 ean 
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Equation (1-11) is the differential equation which 
describes the rate of heat flow through the semiconductor 
Material, with the geometry of Figure 2 and the conditions 
stated before, when a temperature gradient is impressed and 
an electrical current is allowed to flow. 

This equation is non-linear and at least a rigorous 
analytic solution is not obvious. 

If in the range of temperature, when the thermoelectric 
couple is operated, the characteristics of the semiconductor 
material are constant or reasonably can be assumed to be 
constant, then equation (1-11) is greatly simplified and 
takes the following form: 

Mae A dae, Le 0 


+— + = 
de® CT SC Yate? 











eal 


Since the temperature gradient under which the thermo- 
electric element of a cooler is operating is small, equation 
(1-12) reasonably can be used as the basic equation for the 
design. 

This is not the case for the thermoelectric generator 
where the range of temperature is quite extended, and the 
material characteristics present sufficient deviations from 
the corresponding values at low temperatures. 

For this, equation (1-11) will be used as the basis for 


the thermoelectric generator design. 
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An analytic solution of equation (1-11) is given in 
APPENDIX A under certain assumptions. APPENDIX B provides 


the analytic solution of equation (1-12). 


WV, 





III. THERMOELECTRIC COOLER 


A. ANALYSIS OF THE ELEMENTARY COOLER 

The operation of the thermoelectric cooler is based on 
the Peltier effect, which was discovered in the eighteen- 
thirties. This phenomenon was explained correctly by Lenz 
in 1838. When an electric current I flows through a contact 
between two conductors, a certain amount of heat is either 
evolved or absorbed at the contact, depending on the 
directions of the current. This heat is known as the Peltier 
heat Qn: 

Becquerel and others established that Q, is directly 


proportional to the current I: 
CLS uni (3-1) 


Bere (7) is the so-called Peltier coefficient. 

It has also been shown that if the current flowing 
through the junction has the same direction as heat flow 
which is due to the heating of this junction then the 
Peltier heat is absorbed; while in the opposite case, heat 
is evolved. 

The Peltier effect has been observed in a great variety 
of materials, and the Peltier coefficient has been found to 


be dependent on temperature in a complex fashion. 
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The thermocouple is considered under steady state condi- 


tions with its hot junction at a temperature T,, and the cold 


junction at To: 
Electric current is considered flowing through the couple 
in the direction opposite to its Seebeck voltage. 


The arrangement is shown in Figure 3. 


The law of conservation of energy gives: 


2 a te a (3=2) 


where qo is the rate of heat absorption at the cold junction, 
dy is the rate of heat which is evolved at the hot junction, 
and W is the electric power. 

An important efficiency index of the thermocouple is the 


coefficient of performance which is defined by: 





COP = (323) 


Yo 

W 
Another important index is the heat pumping rate — or 

the refrigerating effect =e which can be expressed by: 


1 <TR +, Ans ODA, ee er 
; | 


t=R re Ry 


The thermoelectric coolers are usually operating under 
a small temperature difference. This is particularly true in 


the air conditioning applications. Under these small 


Ee) 








1c) 


FIGURE 3 
ELEMENTARY THERMOELECTRIC COUPLE 





temperature ranges, the variations of the properties of the 
semiconducting materials are limited, and the use of average 
values introduces slight errors of the order of one percent. 

The use of average values simplifies the analysis of the 
system and they are used throughout this section. 

From APPENDIX B the rate of heat flowing to the cold 
junction, because of the existing temperature gradient and 
the generated Joule heat, for the P-type and N-type segments 
as: 

. 5 tal) 
Ft = ~}, Ce poe i : fate) (3-5) 
re Ry tn (2) (+7)p 


Re 
_ & (FT)w A Ge ) : 
TAO Te | = De Ge AT areas, oo 


cB 


Introducing the following variables: 


a(n.) f (3-7) 


a (Ge ip 


tn (33) _ 


as (4 Ty 


%)= d(T) ( la.(7)| = {oC Te)| \ ale (3-9) 


21 





equation (3-4) becomes: 
a % = ar ke = 
Vor a(T)IT. - a + Ae] AT - sl [ x+ey] (3=1G) 


The work done by the electric source, per unit time, 
consists of the Joule heat evolved in the thermocouple arms, 


and the work done against the thermo-e.m.f TAAT 
We TAAT+I°R (3-11) 


where Ris the ohmic resistance of the couple. The coeffi- 


cient of performance can be written as: 


n(rr-( 9+) OT- 31° LBB) 
TAAT +I*(bx+P¥] 


COoP= - (3-12) 


Gore. wa (3-13) 
TAATHIR 


B. CONDITIONS FOR MAXIMUM HEAT PUMPING RATE 

The heat pumping rate of the elementary thermoelectric 
cooler, as it appears in equation (3-10), is a function of 
X, Y, and I. Relative maximum can be found by the Lagrange 


multiplier technique using as constraint function: 


Re 0x + OY (3-14) 
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Thus the following conditions must be satisfied: 


io) 





4+ h® =0 








5 (3-15) 
ie +he=0 (3-16) 
Io _ 

-_"* O (3-17) 


where h is the Lagrange multiplier. 


Eliminating h from equations (3-15) and (3-16) the first 


Condition for °F to be maximum is found to be 
a AY 
i} 
2 | + li 
J) eels 
By substitution of equations (3-7) and (3-8), equation (3-18) 


(3-18) 


becomes: 


(bt) mu _ toe \" 
(PT), Canny 


(Soi) 


The second condition for the maximum heat pumping rate is 
found from equation (3-17) 


ran) 
lee 2 = 
—_ i 





(3220) 


The voltage to be applied to the thermocouple is given by: 


V= o male: Toot (3-21) 
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Substituting the expression for I from equation (3-2) in 


the above equation it is found: 
Vea Tu (3-22) 
The required electrical power for the couple is: 
P= VI = ATy Looe . (3-23) 


The maximum heat pumping rate Ge. is found from equation 


(3-10) after the substitution of equations (3-18) and (3-20). 


Ree Cy Gora) 


Venn” 2R (p-1)(D+t) (3-24) 
where 
tf 
D = |1+3| (3-25) 
g a cme Tutte (3226) 
a ae - 
\/F), 5 Fe5, | 
a eal (e227 
J Te 


24 





The coefficient of performance under the maximum heat pumping 
rate conditions is found from equation (3-13) after trans- 


forming it by means of equations (3-18) and (3-20). 


CoP Se COED: ( DE) (3-28) 
23 (p-1)Cd+4 ) 
The rate of heat rejection is determined with the help of 


equations (3-3), (3-28) and (3-24). 


ele Greece cor +23 (3-29) 


Vs eR (p-4)(D+4) 


C. CONDITIONS FOR MAXIMUM COEFFICIENT OF PERFORMANCE 
As it was in the case of maximum heat pumping rate, the 


coefficient of performance equation (3-13), is also a function 


Gex, Y, and I. 


Applying the same technique, the following conditions 


must be satisfied: 
ok + Le =O (3-30) 
9 x 


—_ 3 he =o (G23) 


OK = (3-32) 
91 
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Where the notation K instead of COP, for the coefficient of 
performance is introduced for convenience, and h again 
represents the Lagrange multiplier. Elimination of h between 
equations (3-30) and (3-31) gives exactly the condition of 
equations (3-18) and (3-19). 

By algebraic manipulation the following expressions are 


obtained: 


(r= a wns (3-33) 
-{ 
_ _AATD_ (3-34) 


(p-1 ) 
TH 


7 te Sar : (3-35) 
© OFmax Bt tL Daa 


It is obvious that equation (3-25) goes to zero when: 


Vu 
r= a 


and it goes to infinity when: 
AT=0 (3-37) 


Combining equations (3-10), (3-18) and (3-33), the heat 


pumping rate becomes: 


a Ty- le Dea) (3-38) 
Vo R(p-1) (De!) . 
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D. COMPARISON OF THE EXTERNAL CONDITIONS 

The two cases of maximum heat pumping rate and maximum 
coefficient of performance differ very markedly as regards 
both the working conditions — voltage required — and the 
efficiency indices. 

Comparing the values of eae and Qo in Figure 4, it is 
seen that ° fae is greater than qo and consequently, in order 
to have the same heat pumping rate for thermocouples working 
under conditions of Gmax and COE aes it is necessary to use a 
larger number of thermoelectric units in the second case than 
in the first. This aspect is not decisive because from the 
efficiency point of view it is first of all necessary to 
ensure that the coefficient of performance is sufficiently 
high, since it is this coefficient which governs the elec- 
trical power consumed in a cooling device. On the other hand, 
it follows from Figure 5 that the coefficient of performance 


at the COP. is greater than COP and that this difference 


ax 
can not be compensated at all. The preferable conditions 
are those when the coefficient of performance is a maximum 
since the materials known at present do not give such a high 
coefficient of performance that it can be afforded to have 
it lower, even for the purpose of economy of material. 
Designs based on the maximum heat pumping rate could 
be used only for the rapid cooling of an enclosure or for 


cooling it to a lower temperature at the cost of additional 


electrical power consumption. 
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E. ESTIMATION OF TEMPERATURES OF OPERATION 
Complete air conditioning involves the simultaneous 
control of the following factors: 
1) Temperature 
2) Humidity 
3) Air motion 
4) Air distribution 
5) Dust 
6) Bacteria 
7) Odors 
8) Toxic gasses 

Of these the first three, temperature, humidity, and 
air motion, are paramount, and no air conditioning system 
deserves the name unless it satisfactorily controls those 
three factors. 

In the present work the thermoelectric air conditioning 
is examined only from the standpoint of the air cooling 
process. 

The device should be also considered as equally appli- 
cable for air heating process by just reversing the direction 
of the current through the thermoelements. 

From the standpoint of human comfort, the object of the 
air cooling device is not to cool the body but to control 
its rate of heat loss within reasonable limits by methods 
that are known to be the most conductive to comfort. 

There are four factors which, by their combined effect, 


determine the cooling power of the atmosphere upon the human 
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body. They are temperature, humidity, air motion and 
radiation. 

The effect of humidity and air movement are matters of 
common experience. When the indoor air is very dry, a 
higher temperature is necessary for comfort because of the 
greater cooling effect of the dry air. The evaporation of 
the perspiration with which the human body is always par- 
tially covered has a marked influence upon the rate of heat 
loss, and the rate of evaporation is dependent upon the 
humidity of the surrounding atneneneees, 

The cooling action of air motion upon the skin increases 
heat removal by the increased evaporative effect and by an 
increase in the convection effect as the layer of air in 
contact with the skin is repeatedly renewed. 

Loss of heat by radiation from the exposed body and 
clothing depends on the temperature of the surrounding walls 
and is independent of air temperature. 

Air movement reduces to some extent the radiation loss 
by lowering the surface temperature of the body. The rate 
of heat loss varies according to Steffen's fourth power law 
and is dependent upon the difference in temperature between 
the surface of the body and the mean temperature of the 
environment. This latter temperature is referred to as the 
mean radiant temperature. 

Experiments made on a large number of human subjects 


have demonstrated that there is a consistent relation between 
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the temperature, humidity and air motion in their effect 
upon comfort. When the results of such observations are 
plotted on a psychometric chart, it is found that the points 
representing conditions of equal warmth fall in approximately 
straight lines. These straight lines are called “effective 
temperature lines". For moving air the effective temperature 
for a given temperature and humidity is lower than that corres- 
ponding'to still air. The effective temperature is the 
only true index of human comfort that is now available. 
Neither the dry-bulb nor the wet-bulb thermometer alone 
trully indicates the conditions of the air from a comfort 
standpoint. 

The optimum conditions for human comfort at summer time 
are considered to be 71 F° effective temperature and air 
movement of 15 to 35 feet per minute for an outdoor temperature 
ef 80 F°. 

The following table gives the corresponding effective 


temperatures for several outdoor dry-bulb temperatures. 


OUTSIDE DRY-BULB EFEECTIVE 
TEMPERATURE (F°) TEMPERATURE (F°) 
105 ee) 

100 75 
95 74 
90 715 
85 72 
80 zal 

TABLE 1 
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These effective temperatures are the statistical results 
of extensive experimental work, where several factors were 
considered, between them the stress which a human body suffers 
when it is moving from low to high temperature environments 
and vice versa. 

The outdoor conditions to be assumed in the calculations 
of the cooling and dehumidifying load should not be the 
highest recorded temperature in a location which usually 
exist for only a few hours during the year. They should be 
chosen as conditions which may be exceeded for a few days 
during the summer. 

In every large city a certain design dry-bulb and wet- 
bulb temperature which are readily ascertained, have become 
established and generally adopted. The following table gives 


these temperatures for a few cities of the U.S.A. 


LOCATION DESIGN TEMPERATURES 

DRY-BULB WET-BULB 
New York City 95 75 
Chicago 95 ID 
Philadelphia 95 78 
Detroit a5 75 
Cleveland 95 75 
se. Louis 95 78 
Washington 95 78 
New Orleans 95 79 
San Francisco 90 65 
Dallas 100 78 
Jacksonville 95 78 

TABLE 2 


33 





Since an automobile is not designed to operate ina 
certain place, the worst case outdoor conditions are selected 
as those corresponding to the city of Dallas. 

So, Outdoor dry-bulb temperature is 100 F° and indoor 


Maintained-effective-temperature the corresponding 75 F°, 
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IV. COOLER DESIGN 


A. GENERAL CONSIDERATIONS 

The design is based on the worst case which corresponds 
to an outdoor temperature of 100 F° and an effective tem- 
perature, for human comfort, of 75 F°. This means that the 
temperature of the incoming air in the hot side of the 
thermo elements would be 100 F° and the maintained tempera- 
ture in the passenger's compartment of the vehicle would 
be 75 F°, 

The proposed configuration of the thermoelectric cooler 
is shown in Figures 6a, 6b. It consists of radially 
arranged thermocouples, forming rings. A number of such 
rings are assembled to form the annular configuration. Fins 
have been introduced in both sides of the thermocouples to 
increase the rate of heat transfer. The outer annular finned 
space forms the heat absorbing exchanger and the inner cylin- 
drical finned space forms the heat rejecting exchanger. The 
cooling medium is air which flows by means of fans, one for 
every heat exchanger. When the vehicle is moving with a 
speed higher than a certain level, there is no need for 
Operation of the heat rejecting exchanger fan. 

It is also assumed that the vehicle has been parked for 
a long period under the sun and the temperature of the indoor 
air has reached 100 F°. 

The temperatures of the hot and cold side of the thermo- 
£ 


elements, T can be estimated if the desired logarithmic 


sl tele 
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FIGURE 6a 


CROSS SECTION OF T.E. COOLER 
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FIGURE 6b 


LONGITUDINAL CROSS SECTION OF T.E. COOLER 
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mean and output temperatures of the cooling air, in both 
sides of the thermocouples are selected. 
According to the definition of mean logarithmic 


temperature: 


= tee - en -~ 
eo oi wes 
i Qanin 


the hot and cold shoe temperatures of the thermocouples are 


estimated as: 


Teg —A tout = Tin (4-2) 
A-4 


where: 
AT 
A =@k 
AT = Tin - area 


K desired value of logarithmic mean temperature 


= temperature differences as shown in 


Bear 0 
ae aS Figure 7. 


The cooling capacity of the thermoelectric cooler is 
selected to be 7000 BTU per hour. This figure is considered 
aS a moderate value, sufficient to serve a small vehicle 
carrying four passengers under the climatological conditions 


described previously. 
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FIGURE 7 


ASSUMED TEMPERATURE DISTRIBUTION FOR 
DETERMINATION OF LOG. MEAN TEMPERATURES 
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As the most suitable thermoelectric material for cooling 
application selected a bismuth — tellurium — antimony com- 
pound for the P-type semiconductor, and a bismuth — tellutium — 
selenium compound for the N-type. These materials have been 
developed since 1968 from Westinghouse. The average proper- 


ties of the materials are displayed in Table 3. 


ITEM UNITS P-type N-type 
Seebeck Coeff. V/deg K 205 x lo ° - 220ieanoN: 
Resistivity ohm-cm 1 28eio* 1.08 x 10°° 
in. COonductiv. watts/deg c-cm 0.0121 0.0149 
Fig. of merit l/deg K Pec Om 3u00 x Cm 

TABLE 3 


In order to obtain good mechanical strength the materials 


are processed by pressure and sintering. 


B. COOLER DESIGN BASED ON MAXIMUM COP 

This design is based on the results obtained from the 
analysis in Section III-c. 

1. Thermoelectric Couple 

The coefficient of performance, as was defined in 

Section III-C, is independent of the geometric configuration 
of the thermoelectric couple. It is only a function of hot 
and cold shoe temperature and of the material properties. 


The heat pumping rate and the optimum current depend on the 
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total electrical resistance of the thermocouple which is 
a strong function of geometry. From this standpoint, an 
optimization of the geometry had been attempted using 


iterative computer techniques and the selected dimensions 


in Table 4. 
ITEM UNITS P-type segment N-type segment 
Ry cm (in.) TO (25756) 7.0 (27756) 
Ro em (in.) 8.5 (3.346) 8.5 (3.346) 
<p degrees 10.0 Unive 
T em (in.) 1.0) .(035937) 1.0 (0.3937) 
TABLE 4 


The calculated cooling and electrical characteristics 


of the thermocouple are shown in Table 5. 


ITEM UNIT VALUE 
Electrical Resistance ohm 2.32 x Lowe 
Hot shoe temp. Th degr. K 322.73 
Cold shoe temp. To degr. K Paps | lama 
Electrical current amperes 14.4 
Applied Voltage volts 0.048 
COP dimensionless Lvisas4 
Required Power watts 0.691 
Heat pumping rate watts 0.74466 
Heat rejection rate watts 1.43569 

TABLE 5 
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For the cooling load of 7000 BTU per hour a total 
number of 2754 thermocouples is required. If a nine degree 
wedge of insulation is provided for each thermocouple, 
then 12 couples form one ring and the total number of rings 
is determined as 230. Each ring, according to Table 4, 
has width 1.0 cm. Allowing insulation of 0.3 cm between 
the rings, the total length of the configuration is 299 cm. 

Four units, each having a length of 74.75 cm and 
consisting of 57.5 rings, will satisfy the requirements. 
The externally supplied electrical power is calculated 


to be: 
2754 x 0.691 = 1903 watts 
and the heat rejection rate: 
2754 x 1.43569 = 3954 watts or 13504 BTU/hour 


2. Cold Side Heat Exchanger 
The introduction of extended surfaces, in the form 
of fins with rectangular cross section, increases the 
exchange surface of the heat exchanger. The fins have an 


efficiency given by: 


z {wh Cmb) (4-3) 


"'Y an 
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where 


m= J2h (4am 


er) 
Figure 8b shows the variation of Np as a function of 
for several values of h. The cooling medium is air with 


average properties displayed in Table 6. 


SYMBOL ITEM UNITS VALUE 
A Thermal conductivity BTU/ft-hour-F° 0.015 
} Dynamic viscosity lbm/ft-sec. Les ome 
p Density lbm/ft? 0.075 
Cp Thermal capacity BTU/lbm-F° 0.24 
TABLE 6 
Figure 9 shows the heat exchanger arrangement. The long 
fins form four different types of air ducts, namely A, B, 
Gy and D. 
Each air duct has an equivalent hydraulic diameter 
given by: 
_ 45 =, 
Doe (4=5) 


where S is the cross section and P the perimeter of the 


euct. 
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FIGURE 8a 


LONGITUDINAL FIN ARRANGEMENT 
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FIGURE 8b 


FIN EPFICIENCY 
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FIGURE 9 


ELEMENTARY SECTION OF COLD SIDE HEAT EXCHANGER 
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The calculated equivalent diameters for each duct 


type are: 

(D.). = 48.284 x 10> ft 
eA i 

(D.}). = 43.89mem108- FE 
e°B és 

(D_). = 43. 60meee Om Ee 
eC F 

(Dey = 62. 00m lon ee 
e’'D rm 


The selected temperatures for the air and thermo- 


couple cold shoe are displayed in Table 7. 


ITEM TEMPERATURE 
Incoming hot air OO Fe 
Outgoing cold air 70 F° 
Maintained room environment TSSES 
Thermocouple cold shoe 2900 Sek 

TABLE 7 


For operation at low noise level, air velocity is selected 
to be 1800 feet per minute. 


The Reynold's number is determined from: 


De V a 
(Nee ED (4-6) 


iF 
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Using the air properties of Table 7 the corresponding 


value of Nor for each type of duct is calculated as: 


(Nop) a = 9285 
(Nop p = 8440 
(Neplc = 8440 
(Nee) p = 11973 


These values are indicative of turbulent flow. The calcu- 
lation of the convection heat transfer coefficients is 


based on the Dittus and Boelter equation: 


0.8 4 0.3 (4-7) 


N = 0.0265 (Nog) pr 


NU 


where Nyy and Nor are the Nusselt's and Prandtl's number 


correspondingly, given by: 
New = h De a 
No. = cet (4-9) 


For a Prandtl's number almost equal to 0.7 the Nusselt's 


number, for the different type of air ducts, is: 
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(Ni) 35.1428 


NU A 
(Niu) B = 32.56 
(Niue = 32.56 
(Niu) p = 43.07 


The corresponding convection heat transfer coefficients are 


found from Equation (4-8) as: 


ha = 10.92 

hh = 11.23 
ha = 11.128 
Ay = 10.376 


Using the above values for the convection heat transfer 
coefficient and assuming aluminum as the fins material, with 
thermal conductivity 92 BTU/ft-hour-F°, the efficiency of 


the long and short fins, according to Equation (4-3), is: 


= 0.9256 and n = 0.9716 


Ney fs 


correspondingly. 
The rate of heat absorption is calculated from 


the following equation: 
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‘= AT 2 hy Sj (4-10) 


where the subscript j corresponds to different types of air 
ducts of the heat exchanger. Analytically the exchange 


effective surfaces per thermocouple column are: 


: 7 2 _ 2 
S, = (Sng, + 3ng, + 1.08) 230 = 1983.2 cm* = 2.135 ft 
2 2 
S, = (Sng, + 0.26) 230 = 1124.24 om* = 1.21 ft 
Pee Say 30 = 1064044 ons a uldde fen 
c al 
Mee (Sn. + Sn. + 1.232) 230 = 201a08eno Olly one 
D fl fs ; 


From the temperatures of Table 2 the logarithmic mean 
temperature of the cooling air, according to Equation (4-1), 


is found to be: 


ATL. = one 


Substitution of the determined values for Ss. and hs in 


Equation (4-10) gives: 


Ca— ceo 20.92 + Loa ierei2zs + 1.146°1b.126 


+ 2.172-10.376) AE, = 1081 BTU/hour 
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For the 12 thermocouple columns of the device, the capacity 


of the heat absorbing exchanger is found: 


Gnor = 12792 BTU/hour 

The above figure for the heat exchanger's capacity offers 

a Margin factor of 1.85 which ensures satisfactory operation 
under even, not considered, worse conditions. 


The air flow requirements can be determined from: 


Se g 
m = — (4-11) 
as Lie vOut 


for ¢ = 12,792 BTU/hour, co = 0.28 BTU/lbm-F°, 


ten = 100 F°, and ore = 70 F°, the rate of air flow is 


determined as: 
m = 1522.8 lbm/hour 
For the air density f = 0.075 lbm/ft® 
: 3 
m= 5.64 ft/sec 
3. Hot Side Heat Exchanger 
This is also a finned type heat exchanger. It con- 
sists of a tube with radially introduced fins of two differ- 


ent lengths. The fin's material is assumed to be aluminum 
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with thermal conductivity 92 BTU/hour-ft-F°. 


Following the 


same procedure as in the cold side heat exchanger the obtained 


results are displayed in Table 8. 


ITEM 


Required capacity 
Hot shoe 
Incoming air 
Temperature 
Outgoing air 
Log. mean difference 
Tube diameter 
Long fin length 
Short fin length 
Number of air ducts 
Duct equivalent diameter 
Duct length 
Air velocity 
Reynold's number 
Nusselt's number 
Prandtl's number 
Convection coefficient 
Long fin efficiency 
Short fin efficiency 
Calculated capacity 
Security factor 


Volumetric rate of flow 


UNIT 


BTU/hour 


Fo 
Fe 
Fe 
F° 
Et 
Et 
Ee 
det 
Et 


ft/sec 


BTU/£t*-hour-F° 


BTU/hour 


£t?/sec 


VALUE 


13,503 
E2des5 
100 
110 

16 
0.4462 
0.1833 
0.082 

24 
0.03784 
955 

86 
20,860 
G75 

0.7 
260 
0.6296 
0.882 
18,983 

1.4 

9.4 





TABLE 8 
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FIGURE 10 


ELEMENTARY T.C,-HOT SIDE 
HEAT EXCHANGER 
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V. THERMOELECTRIC GENERATOR 


A. ANALYSIS OF ELEMENTARY GENERATOR 

The elementary thermoelectric generator, which is operating 
under steady state conditions, consists of a thermocouple 
externally connected to a resistive load, as it is shown in 
Figure ll. 

The operation of the device can be described as follows: 
the hot junction of the thermocouple is assumed at temperature 


Wh The thermocouple arms are assumed to be adiabatically 


c: 
insulated, and heat exchange takes place only at the junctions. 
At the hot junction heat Io is absorbed, and heat qq is 
evolved at the cold junction per unit time. Under these 
conditions, a thermo-e.m.f. E = «AT is established in the 
circuit and it maintains an electrical current I. The 


coefficient in the thermo-e.m.f. is the Seebeck coefficient 


of the couple 





From the law of energy conservation it follows: 


ae ae ‘ ee) 


where W is the available electrical power at the external 


load. 
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FIGURE If 
ELENENTARY T.E. GENURATOR 
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Wie laa 


we TE-I°R (5-3) 
w= TAAT-TFR 


where R is the total electrical resistance of the thermocouple. 
Equations (5-3) suggest the electrical equivalent of the 


thermocouple shown in Figure 12. 


Figure 12: Electrical equivalent of thermoelectric 
couple. 


The general expression for the heat flow rate do has the 


following form: 





‘ye (Ty )1- (0,1) Anst 


n-2,}- POM nals 


T= Tu Vet 





This expression is reduced to a more practical form by 
introducing the estimate of the rate of heat flow in the 


conducting solid, as it is determined in Appendix A 


ae (+ LA ata (% a) a at a 
SAL re] es =} tn(B\ a esse Lé AT (5-5) 


t= Ta 


where 


5 famds 


Lf Ptr) Cr) ( Ly adt) gs 
| [yond | 7 ae 


Tu iu 
4 J 607) ( i AMdT)dt 
a facodt 





e 





OS 


The plus or minus Sign in equation (5-5) allows for two 
possible directions of the electrical current: the minus 
sign representing the current in the P-type segment and the 
plus sign giving the current in the N-type segment. As it 
is shown in Appendix A, the value of e is very close to half 
of the average resistivity for the range of temperature 
Ta! To: 
For convenience in further calculations, it is desirable 


to replace equation (5-5) with a less exact but much simpler 


expression: 




















Alar Rare (PT A 2— a“ = = 
MO TL SA oT a Gia, 16 dT (5-7) 
ToTy = 
where: - 
p>. P(T)dT 


Combining equations (5-4) and (5-7), the rate of heat 


flow at the hot shoe of the thermocouple becomes: 





Sem(ity) let (2D) eo aa 4  & V4, (Pa) _ 
To sili dis 4( >) z (4c), (4tw| z) 
(5-8) 


~1(6,-6,) 7 


Introducing the Kelvin's relation, equation (5-8) becomes: 


got ATT [= Be) aT 2 PL Pex By] TEES) OT Hs 





(yb) 
: ($7) p (4T)y 
and (52100) 
Az (tea (T)-C A, C1)) | AT 
AT Sf N 
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The efficiency of the elementary thermoelectric generator 


is defined as: 


AT AT -1°(@xX+ 20,4] 


= ET Sy) 1a Se ne 
a TyI + [+ xy |ar- $1 [Pre Gl |-1 Cees 


If I and AT are considered as constants, then the efficiency 
N is a function of x and y which determine the geometry of 
the P and N type segments. 

Figure 13 shows the three dimensional representation of 
efficiency as it is interpreted in (5-11), for a certain 


temperature range and S;G, as thermoelectric material. 


B. CONDITIONS FOR MAXIMUM EFFICIENCY 

Equation (5-11) is a function of geometric variables 
xX, y and of the electrical current I, for a given temperature 
range. Using the Lagrange multiplier technique again, the 


necessary conditions for maximum efficiency are defined as: 


9n tee O 

ie = (5-12a) 
Ox = ie 

ot Aine es © ~ (5-12b) 
I ‘ 

cae 6 (5-12c) 


where h is again the Lagrange multiplier and R = 


is the constraint function. 
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FEGURE 13 


THREE DINENSION REPRESENTATION OF EFFICIENCY 
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Elimination of h from equations (5-12a) and (5-12b) gives 


the already found condition 
a ee 
x 
acim res Cw (5-13) 
y Nene 


Equation (5-11) can be written as: 


= 2 
El _alAT -1 R - (5-14) 
X Ty1+ [ar + eV AT -LVR-1(G-H)AT 


Multiplying and dividing by R and introducing the condition 


of (5-13), the efficiency equation becomes: 


= 2 
ay ae 
a ts 
ATwf+ GaT- 4 f- $(6-8y) AT 


where 


i R 
= == —— \2 
= 5-16 
C= (ei, Je5,) (5-16) 
Optimization of (5-15) with respect to f gives the 


condition: 


f a OT 
Gah (5-17) 
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where: 


(5=1e) 
Tav= 4 (T+T) 


The optimum current is: 


«x AT 


R(q+t) (5=19) 


Tort = 


The optimum effective voltage is: 


va AGAT (5-20) 
Goi 
When the conditions (5-13) and (5-17) are satisfied, the 
efficiency reaches its maximum value, which is found from 


equations (5-15) and (5-17): 


= ia @ 
Te 
The electrical power dissipated in the external circuit is 
found from equation (5-19) and (5-20). 


W- XG (AT) 


"Rea al 


The formulae developed above show that the efficiency of 
a thermocouple and its power output increase with increase of 


the working range of temperatures. The thermocouple efficiency, 
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like the cooler C.0.P., depends only the the ratio To/Ty 
and the quality factor G, which is slightly different from 
the quality factor D for the thermoelectric cooler. 

For the present application the working temperatures of 
the thermoelectric generator are varying as the temperature 
of the exhaust gasses varies for different power level 
operation of the automobile engine. In the present state of 
material development, no single thermoelectric material is 
satisfactory for use over the complete temperature range. 
Therefore, it is necessary to use different materials in each 
temperature range. One way to accomplish this, is to form 
each leg of the couple by joining different materials in 
series; that is, forming a sandwich structure of the appropriate 


thermoelectric materials for each temperature range. 


C. OPTIMIZATION OF THE SEGMENTED T.E. GENERATOR 
In Figure 14 two P-type materials, Py and Po, are connected 
in series to form a positive thermoelectric arm and two N-type 
materials, Ny and No, form the negative arm of the thermocouple. 
The optimization problem is to maximize the thermal 
efficiency, subject to specification of the interface tempera- 
tures Ty and Ty over the equivalent length of thermocouple, 


R2 


L = in oe , equivalent cross sectional area Ap=(#T),, 
A 
and the temperature averaged material properties. 


The reciprocal of the thermal efficiency, as defined in 


equation (5-11), is: 


A | RTL + (To-To) Kut Ta-Ty) koa - AIR t Rea) “LEAT (5-23) 
al Taha 
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FIGURE 14 
ELEMENTARY SEGNENTED GENERATOR 
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Introducing the concept of equivalent length and cross section 


for the P-type and N-type segment as: 


c= 
~~ 
~ 
1 
=> 
S 
es 
Paley 
al 


- 


ha - 2) 
buns bn (2) (5-24) 
bare ba 2 
L = Qt L, = Ue 33 Ve 
Ap= (¢T), 
Anz (¢T)n 
The terms in the numerator of equation (5-23) have the 
following interpretations: 
Kes = Aes # 
Mewar ess a (5-25) 
Re B92 
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Th = 
4 J spin) [-9, (dtd 


6. 
Po AT ve 
T Yp,( TdT 
x 
(5-25) 
Ty WF 
23 fut [vided 
ie aes) ee ee 
“AT TH 
bea as ul 
'Y 
The electrical current I in equation (5-23) is given by: 


where: R= Roi + R59 + Ry + Ryo = Rpt Ry is the total 


electrical resistance of the thermocouple and E is the 
circuit voltage produced by the cumulative Seebeck effect 

in the positive and negative elements. The segmented elements 
may be considered as batteries connected in series. 


The voltage E is given by: 


B= gl ae LOE 
where: 
; Ty Ty Tu Ty 
gOS a [4,081 + [ar (T)dT Yo (T)dt -[agoPST| 
* Ve Ty Te 
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ee 


Substituting equation (5-26) into (5-23), the latter becomes: 


E* (Rp, 4 By.) _E_ 
(Ta-Tx) Koy t Ta-Ty) Ky + % Te Bee BBS ~ ate - (n= Ss) “Bae R+Ro (5-28) 


E*R, 
(2+ R,)? 


With the variables 0, and Ge defined by the equations: 


—_ Ty- Tx 
0, — aa, 
lu- Te 
(5-29) 


<7 
By = ie ae 
7 


and multiplying the numerator and denominator of equation 


(5-28) by (R+R,)° , equation (5-28) becomes: 


A. (Ta-T) LOp Key + Ov Ka |CR+R0) + A THE(RtR)-$E*(Rer+Qng) 
E* Re (5230) 


lB 
A A 
_ E CSp4~- Snr) (R+Ro) 


e E* Ro 


With the quantities Rp Ry! Koy: Ky defined by the 


equations: 
/ 
Re= Rede 


Ry = Rw Aw (5-31) 
eee 
Ap 
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equation (5-30) becomes: 


D 


_ (u- Te Ga aoe pee i Tes Ca ils 


N 2 
(5-32) 


4 RQ 
ATy Cae) (Box + Raa) (Ep, Eu) (1+ ay 


ES 
ee ee 


cer ertn= 1.) 2Ro FulGcar. ) 


' | ! ! 
7 (Ty- eye Rot 0, Koa Ry tad a (B2) + Oks Ry [HEN E) 


(5-33) 
- R A oA 
et & Ty (1+ “Re a (Rea Rwva ) - (Sp, ee Gua} (144 =) 
’ ea = 
& (Ta-Te) “ane &' (Tur Te) 
Assuming 7 (4 x | , minimization of equation (5-33) 
with respect _ ot (An) , gives: 
P 
9 (=) > fle ey Po (5-34) 
— Bei Rp > Opihe; Rule) ane 
An An 
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From equation (5-34) the condition is obtained: 


U ) 4 
Bis |Site | 2 (5-35) 
Ae ore Kya Re 


Applying condition (5-35) to equation (5-33), the latter 


becomes: 


i. ee fo Af VEX R 
1 _ [COpReKes) + COw Ry Kwa) 2) C+?) (1+) 
= a a ad 
N E 
(5-36) 
= R a A 

— & Ty Ca) a (Rost Rn) = (6e,- 64) (14-8) 

&' (Ta-T) ZAR A(T) 


With the figure of merit defined as: 


i 3. ' 
es eee a (5-37) 
f 4 ’ ! 
4 f (9, Re Kes) 2 (9, Ri Kn) | 


equation (5-36) becomes: 


42 ote yey ake 


a 4M (ty-7,) 5 oho 1) 


(5-20) 


= (Rp, + Pa) = (6.26 Ge 
2 Ro (Gino te: ) 


So 


With the partial derivative of ea with respect to set 


mn 


equal to zero, the following condition is obtained: 
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4 x 
Tye veo) 


LLL (E04 - Gn - En ) 
qm X' ) 


ye 






(5=29)) 


Substitution of (5-39) into equation (5-38) gives the minimum 


or the maximum efficiency. 


Summarizing, the maximization problem of the segmented 
generator efficiency can be stated as follows: Given the 


Tr gu 


Ss Cc: 
the values of Aye (35. R, Ro: I, and Ls where i= P, N 


values of Kis: eis A As and L determine 


Ti4! 
and j = l, 2. 
Solution of the problem requires that the following 


equations be satisfied: 


! 

% = oo Ci sae) (5-40) 

- RB+Ro 
t t 4 

fe 7 ater re (5-41) 
Ow Kus P 
eed 

: (a) Ty (5240) 


4 | (9p Kea) (Gy 2h by) 2) 
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4 
Perm a) (5-43) 


[= Cos fe Cp = ope bus 


(5-44) 


R= Ror + Rea + Rat Rava (5-45) 
2 2 = = sa 2 y 
41 Rei + <1 Roo -(&,, - %1) Tul - (65,4 = + \pr FE (TT) 5-46) 
ae PC 1x-Te) 


wl a A A ae ae 
: 1 ee + = Rvat Cre Xv) Tal 4 (Ena+ En2)1 ‘ dna ge Ta-Ty) 


= 1G (W- 7) 


(5-47) 


Equations (5-46) and (5-47) represent the energy balance for 
the interfaces between elements Pie Po and Ny. N,- They 
state that the amount of energy flowing to interface, plus 
Or minus the amounts of energy generated or absorbed at the 
interface is equal to the amount of energy flowing from the 


interface. 


D. T.E. GENERATOR DESIGN 

The design of the segmented thermoelectric generator is 
implemented by computer aided iteration methods. Such methods 
have been coded for the IBM-360 as it is displayed in program 


Number one. The procedure consists of several steps: 
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alk) 


2) 
3) 


4) 
5) 


6) 


7) 
8) 


9) 
10) 


For the selected materials, calculation of He e 
Ah) 


rT. : 
am aA ° ° : 
ig Fi for i= P,N and j = 1,2. : 


t 
Calculation of 4,0, Q,, . 


Initial estimation of Gee from: 


pe es a 


— 


pees (7. SN Ae (+ 
Aeon = Doage (eT) 


on oe (Te-Ty) = Ino a (Ty - T) 


NA 


= are bay $ vs 


Calculation of Ry Ry Kpy Kyo , 
oN 

tf 
Ap 


selected value of Ap determination of Ay: 


Calculation of from Equation (5-41) and for 
Determination of figure of merit M from Equation 
(o=37). 

R 


Calculation of = from Equation (5-42). 


Calculation of R from equation: 


R= Ces Qo, ai oa) < (€ Lo x Cxstn2) i 


Calculation of current I from Equation (5-40). 
Calculation of both sides of Equations (5-46) and 
(5-47). If their differences are greater than a 
predetermined figure, then a new set of Log? laa is 
estimated and the procedure is repeated until the 


wanted accuracy is obtained. 
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im) Caiculation of Ray Rye 


12) Calculation of efficiency from Equation (5-38). 


The procedure is based on an "in advance" selection of 


¥ and Ty with the criterion of 


safe behavior of the materials at the extreme operating 


the interface temperatures T 


temperatures. 

When the lengths of the P and N type segments are cal- 
culated for the extreme temperatures, the problem of re- 
estimation of the interface temperatures, for different hot 
junction temperatures, is encountered by modification of 
the iteration procedure in step 10. That is, for the calcu- 
lated lengths, Ty and Ty are varying until Equations (5-46) 
and (5-47) are satisfied. 

The adverse operating conditions of the thermoelectric 
generator — high hot junction temperatures and extended 
operating range — impose restrictions in the selection of the 
best available thermoelectric material for energy conversion. 

All thermoelectric materials represent the best of 
their characteristics in a quite narrow range of temperatures. 
This problem is partially encountered by the segmented type 
generator design, where the one material acts as the comple- 
ment of the other in regions of poor characteristics. 

For the present application, the selected thermoelectric 


Materials and some of their properties are displayed in Table 


9. 
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ITEM MATERIAL MAX OPERAT. TEMPERATURE 








Py Ag Sb Te. 900 (K°) 

Po 25% BioTe, + 75% SboTe, 600 (K°) 

Ny Pb Te (25% Sn Te) 900 (K°) 

No 75% Bi,Te, + 25% Bi5Se 700 (K°) 
TABLE 9 


With proper dopants the maximum operating temperature of Py 
and Ny Materials can be considerably extended. 

The thermoelectric properties of the materials are 
displayed in Figures 15 through 21. The data have been 
taken from Reference [1] and [2] and the corresponding 
temperature relations developed using the "best fit poly- 
nomial" technique as it is shown in Program 2. 

With these data supplied, the computer runs through the 
iteration procedure described earlier, and solves for the 
optimum thermal efficiency and the optimum design parameters. 
The optimum solution is displayed in Table 10. 

The maximum and minimum values for efficiency, current, 
and voltage in Table 10 are associated with the highest and 
lowest hot junction temperature. 

Figures 22 through 24 shows the variation of efficiency, 
voltage, and current with respect to hot Bae od temperature 
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FIGURE 16 
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ITEM UNIT VALUE REMARKS 





Th K? 720 min - 940 max Selected 
To K° 339 Selected 
L= (3) - 0.25845371 Selected 
. Co ln( =) = OF Lo3 
Ti bn (£2) - 0.10545 
Cs = bn (BY) = 0.22275 
Eas bn) - 0.0357 
Ap (eT). rad - cm. a.2 Selected 
AY /AD = 0.5 
Ay= CPT py rad - cm. OL 
R/R - 0.272846 
Ry ohms 1.91 x 107° 
R ohms 7 x Oe 
E volts min 0.135 — maze0-225 
I Amperes min 5.20 9 — maxe23. 6 
W Watts min 0.44 - max 1.41 
7) = min 6% - max 9% 
TABLE 10 
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FIGURE 22 
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E. ESTIMATION OF AVAILABLE HEAT POWER FROM THE EXHAUST GASES 

The recently introduced catalytic converter in the 
exhaust system of automobiles, offers the advantages of 
additional release of heat and higher temperatures, by 
oxidation of the unburned CO and Ho- 

Heat, also, can be extracted by cooling down the gasses 
which come out of the output flange of the converter. 

The thermodynamic analysis which follows is based on 
data provided by FORD motor company. A 2300 om? engine is 
examined and the results of the exhaust gas analysis before 
and after catalytic converter are summarized, for three 
different power levels, in Table ll. 

The analysis is divided into three parts with the neces- 
sary logical sequence. 

The first part contains preliminary calculations, as 
the moles/min and pressure of exhaust gasses. The second is 
involved with the calculation of the available heat energy 
from the catalytic action of the converter, and the third 
is referred to the estimation of the possible extracted 
heat energy by cooling down the gasses after the catalytic 
converter. 

1. Preliminary Calculations 

a) From Table 11 at idle speed, the engine produces 
exhaust gasses at a rate of 10 £t°/min O19 0.2832 m>/min. 
One mole of any gas at the state of 1 atm and 273 K° obtains 


a volume of 0.0224 ai and at standard state — 1 atm, 298 K° — 
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0.02445 m>, Therefore, the flow of the exhaust gasses in 
terms of gram-moles per minute is 11.583. 

Based on the percent mole analysis of Table ll, 
the actual amount of the several constituents of the exhaust 


gas can be summarized as follows: 


CURB IDLE 60 MPH ROAD LOAD PEAK POWER 
GAS gram-mole/min. gram-mole/min. gram-mole/min. 
co, 0.80269 8.2700 MS AUS 
CO 0.4981 0.3475 V2 od 
H50 0500 8.548 24.509 
Ho 0.1865 0.0 5.238 
oy 0.64864 3.5558 729912 
No 8.39707 48.995 198.472 
TOTAL iy. Se3 69.497 266.404 

TABLE 12 


b) The partial pressures of the gasses can be 


calculated from: 


P, = Mx (S240) 


where: 


M, = gram-mole of the gas 
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R = gas constant 
T = gas temperature 


V = volume occupied by all gasses, comprising 
TOT : 
the mixture at temperature T. 


With M, taken from Table 12 the partial pressures 


X 


were calculated for three different power levels, and the 


results are shown in tabular form as follows: 


CURB IDLE 60 MPH ROAD LOAD PEAK POWER 

GAS atm atm atm 

CO. 0.06931 0O.11L90122 0.066007 
CO 0.04301 0.005001 0.048 

H50 0.09067 @.123013 0.09201 

Ho 0.0161 0.0 0.1967 

05 0.056007 0.048005 0.030003 
No 0.725045 0.70508 0.7451 
TOTAL 1.000142 1.0001112 1.00079 

TABLE 13 


Based on the data of Table 13, it is reasonable to 
consider that the chemical process in the catalytic converter 


takes place under constant pressure conditions. 
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2. Released Heat From Converter Action 

The change in enthalpy produced by the procedure 
of a chemical reaction from reactants in their natural state 
at 25° C and 1 atm to products in their natural state at 
25° C and 1 atm is referred to as the "Standard heat of 
meactionatw25° Cc. 

The "standard heat of formation" is the change in 
enthalpy resulting from the formation of a compound at 
25° C and 1 atm from the elements in their natural states 
at 25° C and 1 atm. The "standard heat of combustion” is 
the change in enthalpy resulting from the reaction of a 
compound with elemental oxygen each initially at 25° C and 
l atm to produce specifically defined products at 25° C 
ana 1 atm. 

The heat of a given reaction may be obtained by 
constructing a set of formation and combustion reactions 
such that the sum of the individual reactions is the reaction 
desired. The sum of the enthalpy changes of the individual 
reactions will then be the standard enthalpy change of the 
given reaction. If the heats of formation data are available 
for all compounds involved, the standard heat of reaction 
may be obtained by subtracting the sum of the heats of for- 
mation of the reactants from the sum of the heats of forma- 
tion of the products. Normally, heats of formation are 
available for inorganic compounds and heats of combustion 


are available for organic compounds. 
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Since the initial and final states differ from 25° °¢ 


and 1 atm, the change in enthalpy may be evaluated by a 


series of steps so chosen that the enthalpy change of each 


individual step can determined. 


the following diagram. 


These steps are shown on 


The change in enthalpy for steps 1 through 3 and 


5 through 16 are calculated from: 


AH; = Mi | eid 
4 


where M. is gram-mole and c 
al 
the gas. 
temperature and can be expressed as: 


@_2 
Bape XK +bpT +o, 7 


Maere T is in K°. 


The coefficients a b; and c,; 
Table 14. 

COMPOUND FORMULA Oi 
Carbon dioxide co. 6.214 
Carbon monoxide Co 6.420 
Water H,0 7.296 
Oxygen 05 6.0954 
Hydrogen Ho 6.9469 
Nitrogen No 6.4492 

TABLE 14 


oF 


are displayed 


10.396 
2665 
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Sy Gea eer. 

=O. 1992 
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(5-49) 


is the molal heat capacity of 


The molal heat capacity is a strong function of 


(5-50) 
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CGC. x 107° 
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Step 4 is calculated from tables where the heat of 


formation for different compounds is given. 


The change in 


enthalpy is obtained by subtracting the heat of formation 


of the reactants from that of the products. 


Table 15 displays the change of enthalpy for each 


step, for three power conditions: 


CURB IDLE 
cr Calor/min 
Hy = 1340741 
Ho - 955.36 
H, - 492.60 
Hy, -46303.00 
H_ + 2394.40 
H¢ + 728.60 
Ho_49 + 741.00 
He_13 + 741.30 
Ho_14 + 4859.00 
Hi0-15 + 714.00 
Hi1-16 + 3.40 
TOTAL -38909.37 


60 


TABLE ™ 
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MPH ROAD LOAD 


Calor/min 
- 1418.0 
- 738.7 
00 
-21285.0 
+ 2114.0 
0.0 


930.0 


- 1157.0 


4129.0 


=—S2'62.0 


0.0 


=27.0120-20 


a5 


PEAK POWER 


Calor/min 
- 48025.0 
- 36487.0 
- 20702.0 


-824828.0 


fo 


70449 .0 


+ 


25102.0 


6354.0 


44372.0 


0.0 


+ 885.0 


-892886 .0 





3. Estimation of Available Heat from Hot Gasses 

The gasses flowing out of the catalytic converter 
have a high thermal content, which partially can be extracted, 
introducing a heat exhanger in series with the catalytic 
@eonverter. 

A possible arrangement of converter with heat 
exhanger in series is shown in Figure 25. The change in 
enthalpy, because of the cooling process of the exhaust 
gasses, can be estimated from equation (5-49). 

For a temperature drop 200 F°, the change in enthalpies, 


for different load conditions are shown in Table 16. 


CURB IDLE 60 MPH ROAD LOAD PEAK POWER 
ENTHALPY 
DROP 
BTU/min BTU/min BTU/min 
Eco. 137 10119 Sly 3 
468 
Hoo 5 15 3 
H 102 2400 0 
O 
2 
Hy,0 1355 7846 30700 
He 0 0 1063 
Z 
H 6894 35080 169892 
No 
TOTAL 10090 55460 236296 
TABLE 16 
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CATALYTIC CONVERTER AND HEAT EXCHANGER COMBINATION 
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In summary, the total amount of available heat power 


from the exhaust gasses of the examined engine, under the 


assumed conditions is determined as: 


1) Idle condition 38909 + 10090 50000 cal/min 


= 3486 watts 


2) 60 MPH ROAD LOAD 278264 + 55460 


83286 cal/min 


= 5808 watts 


3) PEAK POWER 892886 + 236296 1129182 cal/min 


78741 watts 


F. FINAL CONSIDERATIONS ON GENERATOR 

The number of required thermoelectric couples for a 
typical design of the generator is based on the available 
thermal power which corresponds to 60 mph road load of the 


automobile engine. 


The design specifications of the thermoelectric generator 


are shown in Table 17. 


96 





ITEM UNIT VALUE REMARKS 


Available thermal power watts 5808 

TH K° 900 SELECTED 
To Ke 35 SELECTED 
Efficiency - 8.9% 

Electrical output power watts S17 

Current amps 22.3 

Voltage volts 25.2 

Number of couples - 43.2 

Rj) Crh. 8.00 

Ro em. Oro 6 

Pp rad 0.2 

oy rad Om 

= em. 1.00 
insulation sais idle 

TABLE 17 


The maximum working temperature of the selected thermo- 
electric materials is restricted to 900 K®°. It is possible 
to have higher temperatures in the wall of the catalytic 
converter when the engine is overloaded. In order to protect 
the hot junction thermoelectric materials from damage, a 
protective system must be provided. Such system could be a 


single by-pass line for the exhaust gasses, automatically 
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operating when a certain temperature is achieved, or intro- 
duction of a heat resisting material between the converter's 
wall and the hot junctions of thermocouples, which will 
insure a sufficient temperature drop. 

The cold junction temperature was arbitrarily selected 
to be 335 K°. It can be even lower, depending on the heat 
rejection system which is employed, i.e. water cooling, forced 
or free air convection. 

The hot junction temperature is varying according to the 
developed power level of the engine. This results in the 
variation of all design specifications of the thermoelectric 
generator, as delivered useful power, voltage, current, etc. 
For this reason, an electrical power conditioning device 
becomes necessary. The function of this device depends on 
the special application of the generator. That is, if the 
generator is employed as the electrical power source of the 
automobile, constant voltage and variable current is required, 
and when the generator supplies part of the electrical power 
to the thermoelectric cooler, constant current and variable 
voltage is desired. 

Contact resistance and conductor losses were not taken 
into account, because the estimation of the former is a matter 
of manufacturing procedure and the latter depends on the 
proper material for the applied contact technique. 

Nevertheless, a reasonable figure of contact resistance 
is of the order of 10 micro-ohms per cm”, which is very small 


compared to the resistivity of the thermoelectric materials. 
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Figure 26 shows the arrangement of catalytic converter with 


the attached thermoelectric generator. 
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FIGURE 26 
T.E. GENERATOR ASSEMBLY 
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VI. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDY 


The obtained figures for the thermoelectric cooler and 
generator could be considered as quite promising at least 
from the standpoint of feasibility of the application. 

The thermoelectric cooler for 7000 BTU/hour was found 
to require 1953.0 watts of electrical power. From this figure 
50 watts are dissipated for the required air blowers. If the 
four units, which comprise the device, are fed in parallel, 
then an external conventional generator, coupled to the 
engine, is needed with 33. volts and 57.6 amperes as rating 
characteristics. 

When the thermoelectric generator and the cooler are 
considered as one system, then the required external electrical 
power is reduced by 18 to 25%, depending on the operating 
power level of the engine, and correspondingly the coefficient 
of performance of the cooler is increasing. 

A direct comparison between the conventional air condi- 
tioning system and the thermoelectric cooler-generator scheme 
is not easy, Siamese there are not available enough data for 
the latter. In Table 18 an attempt for comparison is made in 


the most broad characteristics. 
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CONVENTIONAL THE RMOELECTRIC 


CHARACTERISTIC 

SYSTEM SYSTEM 
Volume Both systems have almost the same 
Weight worst best 
Complexity complex simple 
Maintenance required almost none 
Noise level relatively noisy quiet 
GOP Both systems are in the same level 
Working fluids yes no 

TABLE 18 


From the standpoint of manufacturing cost, the existing 
figures show that the thermoelectric system is two to three 
times more expensive from the corresponding conventional. 

It must be mentioned that the cost of the conventional system 
is referred to the "Series production cost" and the corre- 
sponding for the thermoelectric system, to experimental 
constructions. It is expected that whenever the thermoelectric 
system will be manufactured in series, the cost will decrease 
enough and eventually it will be compatible or even cheaper 
from the conventional system, since the thermoelectric material 
itself is not expensive. 

The analysis showed that the efficiency of the generator 
depends strongly on the thermoelectric material. The obtained 
figures — 6% minimum and 9% maximum — for the generator's 


efficiency, are considered satisfactory compared to 
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efficiencies reported already in the literature. Bates and 
Weinstein developed a similar design using PbTe and SiGe 
which complemented each other over the temperature interval 
from 50C° to 1000 C°. Their generator achieved an efficiency 
of 8.7 percent with a hot junction temperature of 800 C° 
(1073 K°). The present design shows an efficiency of 9% 


with T,, = 940 K° and To = 93 5eK 


H 

The requirements for the primary electrical source of a 
small automobile vehicle are satisfactorily fulfilled from 
the specifications of the segmented type thermoelectric 
generator. 


Table 19 displays the electrical power requirements for 


a small size car. 


ITEM WATTAGE 
(12-Volts) 
high beam 45 
Head lamps 
low beam 40 
Front lamps 21 
Tail lamps 21 
Other lamps 10 
Wiper motor 10 
Electrofan motor ps es 
TOTAL ileeyey 5) 
TABLE 19 
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The thermoelectric generator provides 230 watts in the 
idle condition of the engine and 517 watts in the 60 mph road 
load. So, from this simple examination, the thermoelectric 
generator can replace the conventional alternator, offering 
the great advantages of fuel economy and elimination of 
moving parts. 

The realization of the discussed applications requires 
further theoretical and practical work to be done. The 
systems must be examined thoroughly from the following 
standpoints: 


1) Stress analysis of T.E. couples for T.E. generator 
and T.E. cooler. 


2) Thermal and electrical insulation. 

3) Contact techniques. 

4) Fabrication of thermoelectric pellets. 

5) Power conditioning circuits for T.E. generator. 


6) Control system for the automatic co-operation of cooler, 
T.E. generator and conventional alternator. 


It must be stressed that the present work was developed 
under the scheme of the feasibility study instead of a 


detailed design. 
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APPENDIX A 


APPROXIMATE ANALYTIC SOLUTION OF HEAT BALANCE EQUATION 


The heat balance equation is: 


qonmes ts BOI ane - an sym t Isis it {© - (A-1) 
where: \(T) , @(T), 6(T) . are functions of temperature. 


Equation (A-1) can be written also as: 


ad t T om! _. 
AT) ae a 


ae rah a \oortes * Letrmrst 4 eG} 00a- = 


dT dT ai 
Expression T-— can be written as —— or ——_— 
dr (so) d(dnr) 


If the variable x is introduced such that x = ln r_ then: 


gael. dT 
dr ad x (By 
ag Aa 
r Tay ise 
gfx dpi dry. dg taade . ae 
drldr } > av cant |? Secale ees a 
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2 
ah 4 dt d* 
—— = + 1,21 (A-6) 


de) | Uttdx) Od 
Ae 2 
ue 2. ot = JT (A-7) 
dy? ax? dx 


Applying transformation (A-7) to equation (A-2), it becomes 


ACT) a" [= st. Oo ™ (<t st) + yw é mEscr) St 2 fei: AGpISe 


Or: 


a , Ao) afd eal 
\T) ot (ST aE Neemanmes= . Tact) = 4 
dt dx dx Meas ee) 

Dividing by the product mI equation (A-8) becomes 


me ew (2): CST +12 Lo (A-9) 











The following transformation is introduced 


AM dT (A-10) 


Z=- 
at dx 
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Differentiation with respect to T, equation (A-10) becomes: 


m dlam}] _ m Nap dt dx 


oar | aa @ ‘dxt dt 


or 


aei(sn\_ AG] al am JT (A-11) 


dv ‘dx on dt ae - =e? dx2 


Applying the transformation (A-11) to equation (A-9), the 


latter becomes: 


32 | eee (A-12) 
ze 


dt 


Through the attempted transformation the order of the original 
differential equation was reduced but the non-linearity 
remained, so the solution of it at least is not obvious. From 
this point, the analysis of the problem is based on the 
assumption that in a first approximation the electric effects 
are independent. That means, it is assumed that the Thomson 
effect is not interacting with the Joule effect. 

From this standpoint, the electric effects can be studied 
separately and their influence can be added finally in the 
interpretation of the formulae for the rate of heat flow in 


the two ends of the semiconductor element. 
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1. Influence of Thomson Effect on Heat Conduction 
It is assumed that the Joule effect is not present in 
the system. Then it can be set P(t) = 0. 


Equation (A~-12) becomes, for one possible current direction: 


d2 
AT 





+ Gi) = 0 (A-13) 


Integrating and introducing the transformation (A-10), 


equation (A-13) becomes: 
. -\cry.-m eo tf oenae (A~14) 
where B is the integration constant which must be determined. 


Transforming the independent variable x to the original r, 


equation (A-14) becomes: 
+ 
4 “yar oe Ty) 6Cr) dT (A-15) 
te 


The boundary conditions of the system are: 


Under these conditions, the rate of heat flow in the two ends 


of the element is given by: 
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cot RB I Se a 
~\(T)™ = = 6CT) 
AM “fis Ih TdT 
(A-16) 
4 
_ en —— 
x ) dr l= QR 
Equation (A-15) can also be written as: 
dT B ( ro. d 
Da WG = 1-2 SCIAT ine dy) 
ar — ) (A-17) 
Integrating both sides: 
Ta (BL [1-2 [a-¥ [fscoan] 


The following notations are used: 


Tu 
= ° 34 a0 u 

Q = rate of conduction heat flow = he (2) [ Yond + 

aa Te 

4 

u 

Gp = rate of Thomson heat flow = 1[ seyds 
= 


From equations (A-16) it can be seen that B7Q - 
The equality sign holds for the low temperature end. 
Also, for normal conditions, Q?> | In | and consequently, 


q 
| eal £1, which justifies the required conditions in order 
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to expand the fraction in a series of the 





4-2 [ec 


tOllowing form: 





. 2 T m 
=a age: j2 fecnds\ 4 |< | eee [E fecndt|(a-a9) 
1- % feeds aa oe me 

Br. 


Equation (A-18) now can be written as: 
os Tu + 2 T ~ 
= tT esha = 
ioe Lacnlss( : [senate [ecrdt)-. ( [coxt) Janes 20) 


Integrating term by term and after some algebra, equation 


(A-20) can be written as: 


2 
B=Q{ts ty tay... ty. (A-21) 


where: 


1 _ 1{ zen htt % 
i Goda 


Te 


fom 1 = 232, ..4,n 


From equation (A-20) it can be set, for a first approximation, 


B= Q. Then, B for the second approximation becomes: 
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2 
B= @. 308 Ya... , im (A-22) 





It is sufficient to keep only the first two terms of equation 


(A-22), since higher order terms are quite small. 


Lyf 
=. ACT) eC dtd+ (A-23) 
Ty 
i Ct) dT 
Te 





Ge = Gl. 


Substitution of (A-23) in equation (A-16) yields: 


* 6CT) La@drdt 





“orm rs - 7. > ES) [ 1S 
uz) " ACI T 
(A-24) 
ae yen Lewd d+ 
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re Ro ate [yoo dr 


2. Influence of the Joule Effect on Heat Conduction 


Here it is assumed that the Thomson effect is not 


present, so setting G(T) = 0 equation (A-12) becomes: 


dz ie eEMAG) ae (A-25) 
dy = 
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Rearranging and integrating, it is obtained: 
2 2 Z 
z= D-2 fy e dT 
Te 


where D is the integration constant. 


Introduction of the transformations: 


7 =. AT dT 
3b dx 
ee NS 
= 


in the above equation yields: 





a 
XS DY 


~ (Tmt St =i D, | 4 _~ 20? Lecymer]” (A-26) 
Te 


The appearance of two signs in equation (A-26) leads to the 
conclusion that the temperature distribution along the element 
will present a local maximum Ty, at a point corresponding to 
some radial distance Rat 

Rearranging equation (A~-26), integrating both sides 
and applying the boundary conditions, the following expression 


for D, is obtained: 


(A-27) 


Tam Ty 
pe ) ied \ Ci) cain as + if kee, eS 
INR) aL Benya)’ Z(t 53 Leanne] 


2 





ana for the point (Rie Tn 


= Tan 
Dea: ee eonyend |" (A-28) 


There is a certain value of the current which could 
establish the condition of no back heat flow in the hot end 
of the element. That means, it is possible to have the point 
(Rie Tn shifted to coincide with the point (Ry, Ty) applying 
a certain value of current. 

This value of current can be found from equations 


(A-27) and (A-28) and it is: 


Ty 
a ag, cme AM dT | 
In (2) (A-29) 
tn () ate Se 


The above result might be important for the thermo- 
electric cooler, which should operate under the most efficient 
conditions. 

The rate of heat flow at the two element ends is given 


by the following equations. 








dT 2 arf oe 
Be So es 2 ae OL t- BL emacs] 
(A-30) 
d+ 2) 
pe r=R, 


LES 





It is obvious from equation (A-26), in order for the 


left side to be real, the following condition should exist: 
= a 
D = ne ar oT) AT (A-31) 
] 


for the range of temperature under consideration. 


Also, from equation (A-26) by integration of both sides 


De 2 | * Gian a (A~32) 
2 . a 
GR) La $e [ pcr rda|* 


The condition declared by equation (A~3l1) justifies 
the series expansion of the integrand of equation (A-32), in 


the binomial form. 


T ays 
( ligne af por) \Or) J+) = (rsx) (A-33) 
4° 


is 
2 
where p =~ and x=- aE [ pen yendt 
1 Te 
i 1x2 -1)-2) x3 
Mx) = 44 px+ LCD Eee ee 
Ail) 31 
, PP=4). -- (poman) x™ 
ms 
a x7 45 x3 (an-4) (2n-3)---14 x” 
3 2 eee ic Rts 2p ee ene an 
ee aie ee ae ee 2m mn 
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Application of the above expansion in equation (A-32) and 
integration of it term by term would yield, after algebraic 


manipulation, the following expression for De 


ieee 2 Gye a (-, \* (A-34) 
ta 431 2° (33 ) 


where: 


o 6 abs Jondt 
ae 


n mI 
TH zi 
ep llB) | Lac Ch eoneds Yar ]* 


( [4 oda)" 


Following the same procedure as before, one can set 
D,= Q as a first approximation. Applying this value of D,in 
equation (A-34), D, can be calculated in a second approximation 


as: 


p= Ql t+ eg de 2 w.| (A=35) 


ib Q? 


The higher order correction in equation (A-35) might be 
considered to be quite small to be taken into account in an 


approximation solution. 


1X5 





It is sufficient to approximate D, by: 


D,= aes 


(A-36) 


Applying the above expression for D,in equations (A-30) the 
rate of heat flow at the two ends of the element, now it is 


given by: 


Ty 7 
= yey st | = Bye er ie 2d les) Lac Leourdrdy 
d Yz | Le 
he. (84 e "(faery 


(A-37) 
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Se Composite Solution 


Based on the assumption that the electric effects — 
Joule and Thomson — are non-interacting, superposition of 


the obtained results yield the following composite solution 
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for the rate of flow of conduction heat in the semiconducting 


pellet: 





Ty 7 
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(A-38) 
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(A-39) 


Equation (A-38) is quite complicated to be used in the analysis 
of the thermoelectric cooler or generator, and it can be 


replaced by a simpler but less exact expression as: 














Rp 

. s a s Algesanat Foal) 4g = 
ACT) NT = + [ee | (8 , yor) 5 a Ge (A-40) 

=k, nN te + 

Ra ig 
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a fii —" [sen fyordtet 
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The introduced error in equation (A-40) has been 
investigated using computer techniques and was found to be in 
the order of +5.8% and -3% for the resistivity terms of (A-40) 
for P and N type semiconducting pellets. 

Analytically, a semiconducting SiGe allow has been 
used with electrical resistivity and thermal conductivity as 
shown in Figures 27 and 28. For Ty = 645 K° and To = 400 K° 


the following results were obtained: 


P=2YeE N=TYPE 


5 3 


= 0.64771 x 10. 


COS 
| 


1.12535 x 10. 


“Os 
tl 


3 3 


1 2756s 0m 0.61 x10. 


“Oo! 
it 


i 
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Nie 
“Ol 
tl 


a i 


. /rsitiefeaiee 
e = te te 
( [ acd) 


Te, 


These small errors have a very slight effect in the 
calculated amount of the conduction heat rate of flow as it 
is calculated by (A~38) and (A-40). The calculated error 
was of the order of 0.2%. 

The validity of the assumption of non-interacting 
electrical phenomena, in a first approximation, can be 
justified by replacing the arbitrary relation PO) 1G rere 


with constant values in equations (A-38) and (A-39). Such a 
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FIGURE 27 
THERMAL CONDUCTIVITY x I02 — WATT/CNM/DEG. 
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FIGURE 28 
ELECTRICAL RESISTIVITY - OHM - CM x I02 
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substitution gives the exact form of the rigorous solution 
of the simplified differential equation as obtained in 


Appendix B. 
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APPENDIX B 


SOLUTION OF SIMPLIFIED HEAT BALANCE EQUATION 


The heat balance equation in the thermoelectric element 


with constant properties was defined in Section II as: 


2m 
27; 4 dt B 


ee ee es 


dy? YY dy ce 


same 
Vary 


where B= 


Multiplying by x’, equation (B-1) becomes: 





2 
; 


Introducing the transformation x= l1nr, 


becomes: 


are 
— 


| 


25 


Q 
cal 


Integration of equation (B-3) yields: 


dt = Bx + C, 
Sx 


y= Gy Gx +tBx Sid _ 
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(B-1) 


(B-2) 


equation (B-2) 


(B-3) 


(B-4) 


(B-5) 





Replacing x with ln r equation (B-5) becomes: 
2 
; j2 CQ, +0 nr ~ = 15 (Lar) (B-6) 


The initial conditions for the system are: 


Introducing the above conditions to equation (B-6), it yields: 
2 
Ty = C2 4+Cy Ln Ry 5 * C (den R.) S (B-7) 
2 
Te =C, + Gy hea De -iB Ce 2.) (B-8) 


From equations (B-7) and (B-8) the constants of integration Ci 


and Co are determined as: 


; B) ty (OR as (=)] si 


e 5 —(B-9) 
i (=) 
Se eee & Lee (B-10) 
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The rates of heat flow at the hot junction (T = Ty) and at 


the cold junction are defined as: 


2 fy = = 
bt 4 pe) _. (B-11) 





ES at 
= WiGee ye | = AE 2 & 


= AT — Auli 4 B12 
VCE], TB ney 3 iS a 


The temperature distribution along a SiGe segment is shown 


in Figure 29. 
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FIGURE 29 
TEMPERATURE DISTRIBUTION IN A 51Ge SECNENT 
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